THE LAGRANGE BITOP ON so(4) x so(4) 
AND GEOMETRY OF THE PRYM VARIETIES 

Vladimir Dragovic and Borislav Gajic 

Abstract. A four-dimensional integrable rigid-body system is considered and it is 
shown that it represents two twisted three-dimensional Lagrange tops. A polynomial 
Lax representation, which doesn't fit neither in Dubrovin's nor in Adler-van Moer- 
beke's picture is presented. The algebro-geometric integration procedure is based on 
deep facts from the geometry of the Prym varieties of double coverings of hyperel- 
liptic curves. The correspondence between all such coverings with Prym varieties 
splitted as a sum of two varieties of the same dimension and the integrable hierarchy 
associated to the initial system is established. 
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1. Introduction 



Starting from the end of 60's, the Lax representation has been the one of the 
most powerful tools in the inverse scattering method for the integration of nonlinear 
differential equations, partial and ordinary as well. The Lax equation 



with L{X),A{X) being matrix polynomials in so called spectral parameter A were 
studied in the middle seventies by Dubrovin [12, 15]. That theory was based on 
the notion of the Baker- Akhiezer function, developed by Krichever and others from 
Novikov's school (sec [20, 14]). The L — A pairs of such type and Dubrovin's theory 
were used in the algebro-geometric integration of rigid-body motion in [21, 9]. 

Few years later, Adler and van Moerbeke presented a new version of such a 
theory in [1] . It was based on [22] and also applied in the integration of rigid-body 
motion in several cases [1, 25, 26] 

Recently we have found a Lax representation of that, polynomial, form for the 
classical Hess-Apel'rot system, see [11]. Generalizing it, we constructed a new 
completely integrable system of the classical Euler-Poisson equations of motion of 
a heavy four-dimensional rigid body fixed at a point. Together with generalized 
Lagrange case and generalized symmetric case, which were introduced by Ratiu 
(see [25]), this system exosted the list of integrable Euler-Poisson equations with 
the L operator, which is a quadratic polynomial in A of the form 



The principal aim of this paper is to give algebro-geometric integration of this 
new system. 

However, it appeared that this system did not fit exactly neither in Dubrovin's 
nor in Adler-van Moerbeke's picture. The matrix L satisfies the condition 



Such situation is explicitly excluded by Adler-van Moerbeke (see [1], Theorem 1) 
and not so explicitly by Dubrovin (sec [12], Lemma 5 and Corollary). Up to our 
knowledge, examples which satisfy the condition of type (1) have not been studied 
before. 

Analysis of the spectral curve and the Baker- Akhiezer function shows that dy- 
namics of the system is related to certain Prym variety 11 (which splits according 
to the Mumford-Dalalian theory [24, 10, 28]) and evolution of divisors of some 
meromorphic differentials O^-.Then the condition (1) requires that differentials 



have to be holomorphic during the whole evolution. Compatibility of this require- 
ment with dynamics poseses a strong constraint on the spectral curve: its theta 
divisor should contain some torus. In the case presented here such constraint ap- 
pears to be satisfied according to Mumford's relation (see [24]) 



L{X) = [L{X),A{X)] 



L{X)=X^C + XM + r. 



(1) 



Li2 = L21 = L34 = L43 = 0. 




(2) 



n- c e 
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where 11" is a translation of the Prym variety H. 

The paper is organized as follows. In section 2 the definition of the system of the 
Euler-Poisson equations on so(4) x so(4) is given and few of its basic properties are 
listed such as the L — A pair, a set of first integrals in involution. In the section 3 
transformation of coordinates is performed in classical manner and the connection 
with the Lagrange top is presented. The spectral curve is described in section 
4. In section 5 the Bakcr-Akhiezer function was studied. The next, section 6, 
contains analysis of the Prym variety 11 and via the Mumford-Dalalian theory, the 
connection of algebro-geometric and classical approach from section 3 is established. 
In the section 7 differentials fi* are defined and the holomorphicity condition is 
derived from the condition (1). By using Mumford's relation (2) formulae in the 
theta functions were derived. Necessary gluing of the infinite points of the spectral 
curve and passage to the generalized Jacobian is done in the section 8. The whole 
hierarchy of the Lagrange bitop is considered in the section 9. The higher operators 
Ln, arc polynomials of order N in A. Their spectral curves are double covering of 
hyperelliptic curves of genus 2N — 1 ; these coverings are defined by those divisors 
of order two which are of degree 2N. We can conclude by saying that the Lagrange 
bitop hierarchy realizes all coverings of that kind. 



2. The definition of the system and basic properties 

The equations of motion of a heavy n-dimensional rigid body fixed at a point in 
the moving frame are: 

M = [M,n] + [r,x], 
f = [r,o], 

where the moving frame is such that the matrix / is diagonal in it, diag{I\, . . . , /„). 
Here My = (7^ + G soin) is the kinetic momentum, 17 G soin) is the angular 

velocity, x G so{n) is a given constant matrix (describing a generalized center of the 
mass), F G so(n). Then 7j +Ij are the principal inertia momenta. These equations 

are on the scmydircct product so{n) x so{n) and they were introduced in [25]. 
We are going to consider a four-dimensional case of these equations defined by 

/ xi2 \ 

h=l2=a 000 

h=h = b X34 

V -X34 / 

with the conditions a^h, Xi2,X34 7^ 0, |xi2i 7^ |X34|- 

We will call this system the Lagrange bitop for the reasons we will explain at the 
end of section 3. 

Proposition 1. [11]. The equations of motion (3) under the conditions (4) have 
an L — A pair representation 

|l(A) = [L(A),A(A)] 

L{\) = A^C + AM + F 
A(A) = Ax + 0, 
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where C = {a + b)x- 



Before analysing the spectral properties of the matrices L{X), we will change the 
coordinates in order to diagonalize the matrix C. In this new basis the matrices 
L{X) have the form L{X) = U-^L{X)U, where 





/ 








iV2 
2 


V2 \ 
2 ' 











^/2 


x/2 




u = 




2 




2 








2 


2 




[ 


V2 
2 


is/2 
2 





/ 



After strightforward calculations, we have 



(6) 



where 



L{X) 



/ -iA34 

zA34 

- iPi -i/Ss + Pi 

\iPl+Pl Pl+iPl 



-PI - iPt 
-ipl - PI 

-«Ai2 





Ai2 = A^Ci2 + AMi2+ri2, 

A34 = x^Csi + AM34 + r34, 






Pi = X3 + Xys, 


X3 


(7) 


Pi^ X4 + Xy4, 


X4 




/?3 = S3 + Ay3, 


2/3 




P2=X4 + Xy4, 


2/4 



The spectral polynomial 

p{X, n) = det (L{X) - m • 1) 

has the form 

(8) p{X,^l)=n^ + P{X)n^ + [Q{X)f 



where 
(9) 



P(A) = A^2 + A^^ + Apsp; + 4P4PI 
Q(A) = A12A34 + 2i{p*P4 - PsPl). 



We can rewrite it in terms of and Fy-: 

P(A) = AA^ + BX^ + DX^ +EX + F, 

(10) ^ ,0 

^ ' Q{X) = GX^ + HX^ + IX^ + JX + K. 
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Their coefBcients 



^?i'2 + C|4 = (C+,C+) + (C_,C_), 
2(734M34 + 2(7i2Mi2 = 2 {{C+, M+) + (C_, M_)) , 



B 

D = Mfo + Af: 



23 



Ml'2 



M. 



34 



2Ci2ri2 + 2C34r34 



(11) 



E : 

G 
H 

I : 



(M+,M+) + (M_,M_) + 2 ((c+,r+) + (c_,r_)) , 

2ri2Mi2 + 2ri3Mi3 + 2ri4Mi4 + 2r23M23 + 2r24M24 + 2r34M34 

2((r+,M+) + (r_,M_)), 



34 



(r+,r+) + (r_,r_), 



K = 



12 ^ ^ 13 ^ ^ 14 ^ ^ 23 
C12C34 = (C+,C-), 

C34M12 + C12M34 = {C+, M_) + (C_, M+), 

C34ri2 + r34Ci2 + M12M34 + M23M14 - M13M24 

(c+,r_) + (c_,r+) + (M+,M_), 

M34ri2 + Mi2r34 + Mi4r23 + M23ri4 - ri3M24 - r24Mi3 

(M+,r_) + (M_,r+), 

r34ri2 + r23ri4 — ri3r24 = (r+, r_). 



are integrals of motion of the system (3, 4). We used two vectors M+,M_ e 
which correspond to Mij e so(4) according to 



(12) 



(M+,M_) 





Ml 
-Ml 



\ Ml 



-Ml 


Ml 
M2 



Ml 
-Ml 


M3 



-Ml\ 

y 



Here Ml are the j-th coordinates of the vector M+. The system (3, 4) is Hamil- 
tonian with the Hamiltonian function 



n = -(Mi3f2i3 + Mi4r2l4 + M23f223 + M12O12 + M34f234) + Xl2ri2 + X34r34 

The algebra so(4) x so(4) is 12 dimensional. The general orbits of the coadjoint 
action are 8 dimensional. According to [25], the Casimir functions are coefficients 
of A", A, A"^ in the polynomials [detL{X)Y/^ and -|Tr(I(A))^ 
Since 

1/2 

= GA"* + HX^ + IX^ + JX + K, 



detL{X) 

-^Tr (1(A))' = AX'' + EX + F, 



the Casimir functions are J, K, E, F. Nontrivial integrals of motion are B, D, H, I. 
They are in involution. 

Nontrivial integrals of motion are B, D, H, I are independent in the case Xi2 7^ 
±,)(:s4- When \xv2\ = Ixm]-, then 2H = B or 2H = —B and there are only 3 
independent integrals in involution 

So we have 
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Proposition 2 [11]. For \xi2\ 7^ |X34|> the system (3, 4) is completely integrable 
in the Liouville sense. 

There are two families of integrable Euler-Poisson equations introduced by Ratiu 
in [25] . The generalized symmetric case is defined by the conditions 

/! = ••• = /„, X arbitrary; 

and the generalized Lagrange case which is defined by 

h=l2=a, 73 = ... = /„ = 6, xij=0if(t,i)^{(l,2),(2,l)}. 

The system (3, 4) doesn't fall in any of those families and together with them it 
makes the complete list of systems with the L operator of the form 

L(A) = A^C + AM + r. 

Proposition 3 [11]. If X12 7^ then the Euler-Poisson equations (3) could he 
written in the form (5) (with arbitrary C) if and only if the equations (3) describe 
the generalized symmetric case, the generalized Lagrange case or the Lagrange bitop. 

One can compare this with [24] (Theorem 15, ch. 53) and [16] (Example 2, p. 
1451 ). The Lagrange bitop can be embedded, on the other hand, in the Bolsinov 
construction ([24], Theorem 17, ch. 53). 

3. Clcissical integration and the Lagrange system 

Starting from the well- known decomposition so(4) = so(3) © so(3), let us intro- 
duce 

Ml = ^(M+ -I- M_) M2 = ^(M+ - M_) 

(and similar for Q,r,x), where M+,M_ are defined with (12). Equations (3) 
become 

Ml = 2(Mi X fii + El X Yi) El = 2(Ei x ^1) 
(13) . . 

M2 = 2(M2 x + r2 X X2) r2 = 2(E2 X ^2) 

and 

Xi = (0,0,--(xi2 + X34)) 
X2 = (0,0,-^(xi2 -X34)) 

Since 

M+ = 7+n+ , M_ = 7_f2_ , 
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where /+ = diag(/2 + h, h + ^3, h + h), I- = diag(/i + li, h + h, h + h), the 
connection between Mi, M2 and Vli,Vl2 is 

Ml = + J_)fii + (J+ - 7_)02] 

M2 = ^[(/+ - 7_)0i + (7+ + /_)02]. 

Using (4), we have 

Ml = ((a + 6)0(1)1, (a + 6)0(i)2, (a + 6)^(1)3 + (a - 6)0(2)3) 

M2 = ((a + 6)0(2)1, (a + 6)0(2)2, (a - 6)0(i)3 + (a + 6)0(2)3). 

We denoted with the j component of the vector Oj and similar notation we 
use for r and 

If we denote Oi = (pi,gi,ri), O2 = (p2,92,''2), then the first group of the 
equations (13) becomes 

Pi - mqir2 = -nir(i)2 

(14) qi + mpir2 = nir(i)i 

(a + 6)ri + (a - 6)r2 = 

and 

P2 - mq2n = -n2T(2)2 

(15) 92 + rnp2ri = n2r(2)i 



where 



(a - 6)ri + (o + 6)f2 = 



2(a - b) 2x(i)3 2x(2)3 
m = —, ni = 712 — 



a + a + a + 

The integrals of motion are 

[{a + 6)ri + {a- 6)r2)]x(i)3 = /ii 

(a + bfipl + ql) + [{a + 6)ri + (a - 6)r2]' + 2{a + 6)x(i)3r(i)3 = /12 
(a + 6)pir(i)i + (a + 6)9ir(i)2 + [{a + b)n + {a- 6)r2]r(i)3 = /13 



(16) 



^(1)1 + ^(1)2 + ^(1)3 — 1 

[{a - b)ri + {a + 6)r2)]x(2)3 = /21 
(a - bfipl + qI) + [(a - b)n + (a + 6)^2]' + 2(a + 6)x(2)3r(2)3 = /22 
(a + 6)p2r(2)i + (a + 6)q'2r(2)2 + [{a - b)ri + (a + 6)r2]r(2)3 = /23 

^(2)1 + ^(2)2 + ^(2)3 = 1 

Introducing Pi,ai, defined with pi = picosai, qi = pisinai, using (15) and (16) 
after some calculations, we get 

pj&i + mr2p'i = ni{^^ - air(i)3) 
a + 



(17) 



iipD f = 4nlpl[l - -^(ai + Pin - 4n?(-^ - aiai - ^p?)^ 
nf + ni 



pla2 + mnpl = n2{-^^ - a2r(2)3) 



= ^n\p\{l - ^(a2 + plf\ - 4n^(-^ - a2a2 - '^plf 
ni a + b n2 
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where 

(a + b)ri + (a - b)r2 (a + 6)r2 + (a - h)ri 
"1 = TT "2 = —r 

(a + 

Let us denote ui = , U2 = pi- From (17) we have 

ul = Pi{ui), i = l,2. 



where 
and 



Pi{u) = -Au^-4u'^Bi+4uCi + Di, i = 1,2. 



Bi = 2ai + ai, 

2 2 .O^iX(i)sft-i ,^2, 



(18) a = - a, - - 2a, a., 



Prom the previous relations, wc have 



dui f du2 



So, the integration of the system (13) leads to the functions associated with the 
elliptic curves Ei,E2, where Ei = Ei{ai,ai, X{i)3, fi2, M) are given with: 

(19) Ei : = Pi{u). 

The equations (14) and (15) are very similar to those for the classical Lagrange 
system (see [17]). However, the system (14, 15) doesn't split on two independent 
Lagrangian systems, since the third equations in (14) and (15) together give that 
ri and r2 are constants. Also, in the definition of each of the curves Ei , E2 both 
those constants are involved. That is the reason we refer to the system (3, 4) as 
nonsplitted Lagrange bitop. The formulae (17) are also very close to those for the 
Lagrange system. Although the Lagrange system has a long history (starting from 
1788) and many important papers have been written about it, still there are some 
subtle questions and problems related to its integration (see [16, 5]). So, for the 
further analysis of the Lagrange bitop we pass to the algebro-geometric integration. 



4. The first steps in an algebro-geometric integration procedure 

The L(A) matrix (5) for the Lagrange bitop (3, 4) is a quadratic polynomial in the 
spectral parameter A with matrix coefficients. The general theories describing the 
isospectral deformations for polynomials with matrix coefficients were developed by 
Dubrovin [12, 15] in the middle of 70's and by Adler, van Moerbeke [1] few years 
later. Dubrovin's approach was based on the Baker-Akhiezer function and it was 
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applied in rigid body problems in [21, 9]. The other approach was based on [22] 
and the connection with rigid body problems was given in [1, 25, 26]. 

As it will be shown bellow, non of these two theories can be directly applied in 
our case;. So, wc are going to make certain modifications, and then we will integrate 
the system (3, 4). As usual in the algebro-geometric integration, we consider the 
spectral curve 

r : det (L{X) -fi-l^=0. 

By using (8, 9), we have 
(20) 

r : /i^ + ^2 (a22 + A^^ + 4/33/?3 + 4/34/3:) + [A12A34 + 2i{(3;pi - f^sPDf = 0. 
There is an involution 

a : (A,/i) ^ (A,-/i) 

on the curve F, which corresponds to the skew symmetry of the matrix L{\). Denote 

the factor-curve by Fi = T/a. 

Lemma 1. The curve Fi is a smooth hyperelliptic curve of the genus g(Ti) = 3. 
The arithmetic genus of the curve T is 5a(r) = 9. 

Proof. The curve: 

Fi : w2 + p(a)m+[Q(A)]2 =0, 
is hyperelliptic, and its equation in the canonical form is: 

(21) uj = iM! _ [g(A)]2, 

where ui = u + P{\)/2. Since ^^^^^^ [(^(A)]^ is a polynomial of the degree 8, the 

genus of the curve Fi is giTi) = 3. The curve F is a double covering of Fi, and the 
ramification divisor is of the degree 8. According to the Riemann-Hurwitz formula, 
the arithmetic genus of F is ga(^) = 9. 

Lemma 2. The spectral curve F has four ordinary double points Si,i = 1, ... ,4. 
The genus of its normalization F is five. 

Proof. From the equations 

dp{X, m) ^ Q dpjX, m) ^ Q 
dX ' djj. ' 

the four ordinary double points arc Sk = (A/j, 0), fc = 1, .., 4, where Xk arc zeroes of 
the polynomial (5(A). Thus, the genus of the curve F is .9(F) = (?a(r) —4 = 5. 

Lemma 3. The singular points Si of the curve F are fixed points of the involution 
a. The involution a exchanges the two branches ofT at 5,. 

Proof. Fixed points of the a are defined with ^ = 0, thus 5*^ are fixed points. Since 
their projections on Fi are smooth points, the involution a exchanges the branches 
of F, which are given by the equation 

2_ -P(A) + VP2(A)_4Q2(A) 
^ ~ 2 

In general, whenever the matrix L{X) is antisymmetric, the spectral curve is 
reducible in odd-dimensional case and singular in even-dimensional case. 

Before starting the study of the analytic properties of the Baker- Akhiezer func- 
tion, let us give the formulae for (nonnormalized) eigen-vectors of the matrix L{X). 
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Lemma 4. If the vector f = (/i, /2, fa, f^Y is given by 



A = (A?2 + M2)(iA34 - /i) - 2/x(/J3/?3* + + 2Al2(/?3/?4* ' f^if^t) 

(22) /2 = 2M/33-*/34)(^/33*+/3:) 

/4 = (*/33* + PX) [(«Ai2 + /x)(«A34 - m) + 2^(/33/3: - /33/34)] 
i/ien L{\)f ~ fjtf. 

Corollary 1. The eigenvectors f normalized by the condition 

f[ + f2 + fa + f'i = h 
have different values in the points S[, S'/ G T, which cover the singular points Si G T. 



5. The Baker- Akhiezer function 

The general integration technique based on the Baker- Akhiezer function was 
developed by Krichever (see [20], [14] and bibliography therein). The application on 
the matrix polynomials was done, as we said, by Dubrovin (see [12, 15]). Following 
those ideas, we consider the next eigen-problem 



|-+i(A))Vfe = 0, 



(23) \dt 

L{X)tpk = IJ'ki'k, 

where tpk are the eigenvectors with the eigenvalue fik- Then ipk{t,X) form 4x4 
matrix with components ipKt, A). Denote by (f^ corresponding inverse matrix. Let 
us introduce 

mt,T, (A,/Xfe))=^Ui,A)-^,^(r,A) 

and 

ffj(i, (A, Mfc)) = t, (A, Mfc)) = Mt, A) • ^']{t, A) 
(there is no summation on k) or, in other words 

g{t)=Mt)'^V>{t)''. 
It is easy to check that the function ^j(f, r, {X, Hk)) satisfies 

|- + i(A)) *(t,T, (A,M,)) = 



Then, if we denote with $(t. A) the fundamental solution of 

^+i(A)) $(t,A) = 0, 
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normalized with $(r) = 1, there is a relation 

T, (A, Hk)) = ^t, X)g{T, (A, Hk)). 
Matrix g is of rank 1, and we have 

dip A I df 4 ^5 r 7i 

We can consider vector-functions ipk{t^^) = -^)) V'fc(^i '^))^ as one vector- 

function ip{t, (A,/Lt)) = (ip^it, (A,/u)), (A,/u)))'^ on the curve F defined with 

ip'^{t, (A, /ife)) = ipl{t,X). The same wc have for the matrix iy9*^. The relations for 
the divisors of zeroes and poles of the functions i cpi in the affine part of the 
curve r are: 



(24) {g'j))^ = dj{t) + d\t)-Dr-D',, 

where Dr is the ramification divisor over A plane(see [12]) and Ds is divisor of 

singular points, D'^ < Dg. One can easily calculate deg Dr = 16, degDg = 8. 

The matrix elements g]{t,{\, fjtk)) are meromorphic functions on the curve T. 
We need their asymptotics in the neighborhoods of the points Pfc, which cover the 
point \ = oo. Let ?/'*; be the eigenvector of the matrix L{X) normalized in Pj, by 
the condition = 1, and let (ff be the inverse matrix for tpl. We will use the 
following 



Lemma 5. The matrix elements of g have another decomposition 

g] = i^iif^ = i>i<fi']. 



The proof of the Lemma is an immediate consequence of the proportionality of 
the vectors ipk and ipk (and cp'^ and cp'^). 

Lemma 6. (a) The matrix g has a representation 

_ IJ,^ + aifj? + a2H + as 
^ dp{\ii)/dn 

where 

ai = L,a2 = PE + L^, as = Pi + L^- 
(b) For the Lax matrix L from (3) it holds 

as = P{Xi)L{Xi) + L^{Xi) = 0, 

for Xi : Q{Xi) = 0. 

The proof of the Lemma follows from [12] and strightforward calculation. From 
the (a) part one can see that the matrix g could have poles at the singular points 
of the spectral curve. But, from (b) we have 
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CoroUeiry 2. The matrix g doesn't have poles at the singular points of the curve 

r. 



So, from now on, taking Corollaries 1 and 2 into account, we will consider all the 
functions in this section as functions on the normalization T of the spectral curve 

r. 

Since the functions ipl i ipj are meromorphic functions in the neighborhood of 
the points Pk, their asymptotics can be calculated by expanding •i^/; as a power 
series in j in the neighborhood of the point A = oo around the vector e^, where 



"k — "k 



SI- We get 



C + 



M r 

A2 



Ui Vi Wi 

A2 + A3 



r - -I- — -I- — -I- — 
" + A + A2 + A3 



A +A^+A^ 



where the matrices C, M and f are defined by 

L(A) = + AM + f . 
Comparing the same powers of ^ , we get the system of the equations 



(25) 



Csi — CiiCi 

Cui + Me,; = CuUi + hei 

Cvi + Mui + f Ci = CiiVi + biUi + diCi 

Cwi + Mvi + fui = CiiWi + biVi + diUi + /ije,. 



Prom the system (25), we have 



(26) 



{Ui)i 
{Ui)j 



0, ivi)i = 0, {wi)i 







Mi 



Cii — 


^33 


1 




Cii — 


^33 


1 





E 



3 + i 

MjkMki 



\kjti 



Cii — Ckk Cii — Cj 



33 



Ci — Cj 



^Mjkivi)k + ^i'jk{ui)k - bi{vi)j - di{ui)j 



bi = M„ 



X - MikMki 



hi = ^Mik{vi)k + y^^tjkjui) 

k^i k^i 
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So the matrix ip = {V'fe} in the neighborhood of A = oo has the form: 

Let the expansion of the matrix (p = {(f^} he 

Ul Vi Wi ^ / 1 



From tp ■ (fi = 1 we obtain 



Ml = —u 



(28) vi=u'^ -V 

wi = 2uv — —w. 

Prom g = tph <Si <p'' iJi the neighborhood of A = oo it holds 

(29) g]{tiK^^k))^i>lit,X)^■it,X)■ 
The last relation and (26) and (28) imply 

9j{t, (A, /ifc)) = 1 for i = j = A;. 
For i = j ^ k and i, fc 6 {1, 2} oi i,k G {3, 4} we have 

Qjit, (A,/ife)) = ^Vik{vi)ki + (^-^ 

So the functions have a zero of the fourth order at the point Pk- Here we used 
the fact that M12 = 0, M34 = 0. 

On the other hand, iov i = j ^ k and i € {1,2}, k e {3,4} or fc € {1,2} and 
i G {3, 4} we have 

Thus, the functions g] in this case have a zero of the second order at Pk- 

In the same manner, for i ^ j, k ^ i, k ^ j, i,j G {1,2} or i,j G {3,4}, it holds 



9i{ j{ jfJ- )) ^2 ^Q^^ _ Cii){Ckk - Cjj) \ A^ 

So gj have a zero of the second order at Pk- 

For i ^ j, k = i, i,j S {1,2} or i, j e {3, 4} we have 



fl}(i,(A,Mfc)) = ^ + 0. ^3 
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and have a zero of the second degree at Pk- 

In the case i ^ k = j, i,j £ {1, 2} or i,j € {3,4} the functions Qj at Pk also 
have a zero of the second degree. 

Fori^j,ky^i,k^j and ^ {(1,2), (2,1), (3,4), (4,3)} we have 

7/, /^ ^^ const „ / 1 

g}{t, (A,Mfc)) = ^+0(^^ 

So, the degree of a zero at Pk of the functions gf in this case is 3. 
In the last two cases the functions g^ have simple zeroes at Pk'. 
for i^j,k = j, ^ {(1, 2), (2, 1), (3, 4), (4, 3)} the expansion is 

(30) 5i(t>(A,Mi)) = ^pr^^ + o(^ 

and similarly for j, k = i, ^ {(1,2), (2, 1), (3,4), (4,3)} we have 

1 Mji „ / 1 



Summarizing, we get 

Lemma 7. The functions have the following divisors in infinity: 
fori^j, i {(1,2), (2,1), (3,4), (4,3)} 

{9})x=oo = 3 (Pi + P2 + Ps + Pa) - 2Pi - 2P,-; 

forij e {1,2} ori,j e {3,4} 

{g'j)x=oo = 2 (Pi + P2 + Ps + P4) for i ^ i, 

(s;)a=oo = 2 (Pi + P2 + P3 + P4) - 2Pi + 2Pj for i = j, 
where I ^ i and l,i & {1, 2} or l,i € {3, 4}. 

Let us denote by dj and by rf' the following divisors: 

di=dl+P2, d2 = d2 + Pl, d3 = d3 + P4,, d4, = d4 + P3, 

d^ = d^ + P2, d^ = d^ + Pi, d^ = d^ + Pi, d'^ = d'^ + P3. 

Finally, wc have 
Proposition 4. 

a) The divisors of matrix elements of g are 

(31) (5]) = d' + dj-Dr + 2 (Pi + P2 + P3 + Pi) -P^- Pj 

b) The divisors di,d^ are of the same degree 

deg di = deg d^ = 5. 
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Proof, a) is a consequence of the previous lemma. 

From degdi = deg = degdi + 1 = degd^ + 1 and degDr = 16 using (16) we 
have b). 

For the functions 



mt,T, (A,Mfe)) 



Y.s9j{T-, (A,Mfe)) 
we have 

(32) r{t,r,{X,tik))=J2^i{t,XW{T,{X,iik)) 

S 

where are the eigenvector of L(A) normalized by the condition 

J^h'{t,{X,iik)) = l 

s 

From (32) it follows that 



(33) Vi'(t,T,(A,/Xfe)) = ^$U*,A) 



Proposition 5. The functions ip"^ satisfy the following properties 

a) In the affine part ofT the function has 4 time dependent zeroes which belong 
to the divisor d^{t) defined by formula (24), and 8 time independent poles, e.g. 

{i>\t,T, (A,Mfe)))^ = d\t) - -D, degP = 8. 

b) At the points Pk, the functions ip''' have essential singularities as follows: 

i>\t, T, (A, /i)) = exp [-{t - T)Rk] a\t, r, (A, fi)) 
where Rk are given with 

Ri = i (^+^34) ,i?2 = -i?i,i?3 = i (^+^12) ,i?4 = -R3 

and d' are holomorphic in a neighborhood of Pk, 

ff{T,T, {X,ll))=h\T, {X,H)) 

and 

a\t,T,Pk)=6^ + vlit)z + 0{z^), 

with 

i Mki 

^ii ~ ^kk 
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Proof. From (33) wc see that tp"^ has d^{t) as a divisor of zeroes in the afhne part. 
Also from (32) it follows that poles of tp"^ are poles of /i", and do not depend on 
time. The functions /i* are meromorphic on F, and they have the same number of 
poles and zeroes. 

From (25) and (26), we have that 

-h^ has simple zeroes at P3 and P4, and the double zero at P2; 
-h^ has simple zeroes at P3 and P4, and the double zero at Pi; 
-h^ has simple zeroes at Pi and P2, and the double zero at P4; 
-/i^ has simple zeroes at Pi and P2, and the double zero at P3. 
As in [20, 14], it could be proved that the functions have divisors of poles 
I), deg l) = 8. This proves the rest of a). 
From 

dt Ht,{\,,i)) rkit,x) 

and asymptotics (26), (27) for t/)^ in a neighborhood of the points Pfe, using the 
proportionality of tjj'' and ■ip'', we obtain 

r{t,T)=exp[(t-T)Rk]a\t,T). 

Starting from the expansion in a neighborhood of Rk 

r^t, A) = e«'=*(4 + viit)z + 0{z^)), v'lit, A) = e-«'=*(6*= + wf{t)z + Oiz^)), 

using (30) we get — 51 and 

i Mki 
vi = — —. 

This finish the proof of the proposition. 

Lemma 8. The following relation takes place on the Jacobian Jac(r).- 

Aid\t) + ad\t)) = Aid^T) + ad^T)) 
where A is the Abel map from the curve V to Jac(r). 

Proof. Let us introduce functions (p^{t,T, {X, 11)) = ^i(^'o'(^'ff) ■ Then (y'*)a = 
d^(t) — d^ir). Using the relation 

a(Pi) = P2, a(P3) = P4, 

and statement b) in the Proposition 5 wc see that aip^ -Lp"^ arc meromorphic functions 
(they do not have the essential singularities in P^). Consequently, for their divisors 
of zeroes and poles, it holds: 

{oip' ■ ip'{t,T, (A, At))) = d^{t) + crd'{t) - d^{T) - ad^{T). 

Applying the Abel theorem, we finish the proof. 
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From the previous Lemma we see that the vectors belong to some trans- 

lation of the Prym variety 11 ~ Prt/m(r|ri). More details concerning the Prym 
varieties one can find in [24, 23, 19, 7, 29, 28, 3, 6]. The natural question arises to 
compare two twodimensional tori 11 and Ei x E2, where the elliptic curves Ei are 
defined in (19). 



6. Geometry of the Prym variety 11 

Together with the curve Fi, one can consider curves Ci and C2 defined by the 
equations 

(34) c,:v^ = s{^^+Q{X)) 



(35) C2:v' = - Q{X)) 

where s is a constant to be fixed in the next Lemma. 



Lemma 9. If s = 2/{a + 6)^ in (34, 35) then the curves Ci are such that 

Ei = Jac{Ci) i = l,2. 

Proof. The curves Ei from (19) can be represented in a canonical form 

(36) 2/^ = - 52iM - ff3i, « = 1,2, 

where 

^ 3 ^ ^J' y^27 3 4 ' 

On the other hand, the equations for curves C, are 

y'^ = sia^X^ + afX^ + afx'^ + afX + af), 

where ^ 

a± = (Ci2±C34)(Afl2±Ai34), 

at = ((M12 ± M34)' + (M23 ± M14)' + (Mi3 T Af24)')/2 

+ (Cl2±C34)(ri2±r34), 
at = (M12 ± M34)(ri2 ± r34) + (M23 ± Mi4)(r23 ± ri4) 
+ (Mi3TM24)(ri3±r24), 

at = 2. 

The sign + corresponds to the curve Ci and — to C2. 
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Using the fact that the Jacobian of the curve given with the equation 
2/2 = + 45^3 + 6cX^ + AdX + e 
is a canonical curve of the form (36) with 

g2 = ae — ibd + Sc^, = ace + 2bcd — ad^ — eb^ — c^, 

we finish the proof of the Lemma by strightforward calculation. 

Since the curve Fi is hypercUiptic, in a study of the Prym variety 11 the Mumford- 
Dalalian theory can be applied (see [28, 24, 10]). Thus, using the previous Lemma, 

wc come to 

Theorem 1. The following relations take place: 

a The Prymian U is isomorphic to the product of the curves Ei : 

n = Jac{Ci) X Jac{C2). 

b The curve T is the desingularization of Fi x (P^)C2 = Ci x (Pi)ri. 
c The canonical polarization divisor EofTl satisfies 

H = £i x 02 + Oi x E2 

where 9, is the theta - divisor of Ei . 

Proof. The Prym variety 11 corresponds to the unramificd double covering tt : F ^ 
Fi. This covering is determined by the divisor D G Jac2(Fi), such that 2D = (/x). 
So 

D = Ri+R2 + R3 + Ri- 2(Pi + P3) 

where Ri = 7r(S'i) arc the projections of the singular points on F and Pj = 7r(Pi) 
are the projections of the infinite points on F. 

On the other hand, the double covering over Fi defined by the curves Ci , C2 
corresponds to the divisor Di G Jac2{Ti) 

Di = Xi + X2 + X3 + X4 - 2(Pi + P3) 

where Xi are the branch points on Fi defined by P/2 — Q = 0. Simple calculation 
shows that 

/x + P/2 2 ^ P/2-Q 
^ l^-Q ' -2/x 

holds on Fi. Prom the last relation it follows that the divisors D and Di are 
equivalent. 

The rest of the theorem now follows from the Mumford - Dalalian theory [10, 
24, 28]. 

Theorem 1 explains the connection between the curves Ei , E2 and the Prym 

variety H. Further analysis of properties of the Prym varieties necessary for the 
understanding of dynamics of the Lagrange bitop will be done in the next section. 
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7. Isoholomorphisity condition, Mumford's relation and 
integration using the Baker-Akhiezer function 

According to the Proposition 5, the Baker - Akhiezer function ^ satisfies usual 
conditions of normalized (n=)4-point function on the curve of genus g = 5 with 
the divisor "D of degree dcgV = g + n — 1 = 8, sec [14, 13]. And by the general 
theory, it should determine all dynamics uniquely. The basic question is why is such 
dynamics compatible with the condition (1)? In other words, why is the evolution 
of divisors di{t) such that all the time d\ contains P2, ^2 contains Pi and so on. To 
answer this question, let us consider the differentials OJ- 

VL]=g^d\ i,j = l,...,4. 

It was proven by Dubrovin in the case of general position, that fi* is a meromorphic 

differential having poles at Pi and Pj, with residuums and —vl respectively. 
We have a simple 

Lemma 10. The condition (1) is equivalent to 

(37) vl =vl=vl=vt=0. 

From the Lemma 10 and Corollary 2, it follows 
Proposition 6. The four differentials 

are holomorphic during the whole evolution. 



We can say that the condition (1) (together with the Corollary 2) implies iso- 
holomorphicity. Let us recall the general formulae for v from [13]. 

_ Xi0{A{P^)-A{P,)+tU + zo) 
^ ' '~ X,0{tU + zo)e{P,,P,) 



where U = ^ x^^^ U'^^^ is certain linear combination of b periods f/^*^ of the differ- 
entials of the 
scalar s, and 



entials of the second kinde fip?, which have pole of order two at Pi; Aj are nonzero 



e(P.P,)^ 



-a^(„0[9](o))V2(_a^,,,0[9](o))i/2)- 

(Here 9 is an arbitrary odd nondegenerate characteristics.) 



Proposition 7. Holomorphicity of some of the differentials 12* implies that the 
theta divisor of the spectral curve contains some torus. 



In a case of spectral curve which is a double unramified covering 
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with S'(ri) = g, g'(r) = 2g — l, as we have here, it is really satisfied that the theta 
divisor contains a torus, see [24]. More precisely, following [24], let us denote by 
n~ the set 

n- = |l e Pic^3-'^t\NmL = Kri,h^{L) is odd} , 

where Kp^ is the canonical class of the curve Fi and Nm : PicT — > PicTi is the 

norm map, sec [24, 28] for details. For us, it is crucial that 11" is a translate of the 
Prym variety 11 and that Mumford's relation ([24]) holds 

(39) n- c Op. 
Let us denote 

(40) U = lixuU^'^ - X3iU''^^ + xi2^7(3) _ ^^2C/(4)), 

where C/'^*^ is the vector of b periods of the differential of the second kinde flp\ 
which is normalized by the condition that a periods are zero. We suppose here 
that the cycles a, b on the curve F and a, b on Fi are chosen to correspond to the 
involution a and the projection tt, see [2, 28]: 

(T(ofe) = Ofe+2, A; = 1,2; 

(41) 

7r(ao) = ao; 7r(6o) = 2bo. 

The basis of normalized holomorphic differentials [uq, ■ ■ ■ , U5] on F and [vo,vi,V2] 
on Fi are chosen such that 

(42) -*{vo)=uo, 

IT* (vi) = Vi + a{vi) = Vi + Vi+2, i = l,2. 

Now we have 



Theorem 2. // the vector zq in (38) corresponds to the translation of the Prym 
variety 11 to 11", and the vector U is defined by (4O) than the conditions (37) are 
satisfied. 

Proof. Proof follows from the relations (38) and (40) and fact that P2 = cr{Pi) and 
Pi = aiPs). 

Proposition 8. The explicit formula for zq is 

Zo = 2(''^oo,'roi,'ro2,Toi,fo2), 

where 

Toi = Ui, i = 0, 1, 2. 

Jbo 

The proof follows from [19], the Proposition 4.7. 

The formulae for scalars Aj from the formula (38) will be given in the next 
section. 
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8. The evolution on the generaUzed Jacobian 

The evolution on the Jacobian of the spectral curve, as we considered Jac(r) in 
the Section 7, gives the possibility to reconstruct the evolution of the Lax matrix 
L{X) only up to the conjugation by diagonal matrices. That was one of the limita- 
tions of the Adler-van Moerbeke approach (see [1], Theorem 1). To overcome this 
problem, we are going to consider, following Dubrovin, the generalized Jacobian, 
obtained by gluing together the infinite points; in the present case Pi,P2,P3,P4 : 

Jac(f|{Pi,P2,P3,-P4}). 

It can be understood as a set of classes of relative equivalence among the divisors on 
r of certain degree. Two divisors of the same degree Di and D2 are called equivalent 
relative to the points Pi,P2, P3, P4, if there exists a function / meromorphic on F 
such that (/) = Di-D2 and /(Pi) = /(P2) = /(P3) = /(P4). 
The generalized Abel map is defined with 

i(P) = (A(P),Ai(P),...,A4(P)), 

where A{P) is the standard Abel map and 

Xi{P)=exp Q.p.Q„,i = l,...,A. 

JPo 

Here f^PjQo denotes the normalized differential of the third kinde, with poles at Pi 

and at arbitrary fixed point Qo- 

Then the generalized Abel theorem (see [19]) can be formulated as 

Lemma 11 (the generalized Abel theorem). The divisors Di and D2 are 

equivalent relative to the points Pi, P2, P3, P4. if and only if there exist integervalued 
vectors N, M such that 

A{Di) = A{D2) + 2'kN + BM, 

Xj{Dx) = cXj{D2)exp{M,A{D2)),j = 1, ...,4 

where c is some constant and B is the period matrix of the curve T. 

The generalized Jacobi inverse problem can be formulated as the question of 
finding, for given z, points Qi, . . . , Qs such that 

8 4 

^A(Q,)-^A(P,)=z + i^, 
1 2 

Xj = cexp22 / ^PjQo + «i> J = 1> ■•■4, 

where the constants kj depend on the curve F, the points Pi,P2,P3,P4 and the 
choice of local parameters around them. 

We will denote by Qs the points which belong to the divisor V from the Propo- 
sition 5, and by E the Pryme- form from [19]. Then we have 
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Proposition 9. The scalars Xj from the formula (38) are given with 

Xj = X°exp^ix^'''>-fj, 

where 

8 cQs 



X° = cexpV / ^p,Qo + i^j, 

„_i JPn 



vector X = {x^^\. . • ,a;W) is i(x34, -X34,Xi2, -X12) 

To give the formulae for the Baker-Akhiezer function, we need some notations. 

Let 

where ^ 

JPo 



and 7™ is defined by the expansion 

nP 
'Po 

Denote 

fP 

'Po — ■■■ ^\^y^) -^y^j 

Now wc can state 

Proposition 10. The Baker-Akhiezer function is given by 

yo e(AiP-p,)-zo) 

^{P,p,) 

Z^fe=l '^k €{P,Pk) 



r n^^l = -km + 7Z + 0{k-'), P ^ Pm- 
JPo 



i^^ix^p) ^ <p^{x,p) _^ ' ,„;/r;p-p„.-.„. . j = i,-.-,4. 



The proofs of the statements in this Section are standard from Dubrovin's ap- 
proach. 

Having established how parameters of the formula (38) evolve, the reconstruc- 
tion of the evolution of the phase space variables follows immediately from the 
Proposition 5. The generalized Liouvillc tori arc four dimensional. Since two of 
the integrals of the motion of the Lagrange bitop are linear (see (11)), according 
to the well known fact of Classical Mechanics ([31, 3]) those generalized tori have 
twodimensional afRne part. The two-dimensional compact part of such a torus 
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corresponds to the real part of the two-dimensional Prymian H. The affine part 
corresponds to the odd part of the afRne part of the generalized Jacobian 

Jac(f I {PuP2,P3,Pi}) = Jac(f) X C*\ 

From the Theorem Ic, it follows that the reduction of the formulae can be done up 
to the elliptic theta functions on Ei and exponential functions. 



9. The Lagrange bitop hierarchy and 
equallysplitted double hyperelliptic coverings 

According to the Mumford - Dalalian theory (see [10, 24, 28]), double unramifide 
coverings over a hyperelliptic curve = P2g+2ix) of genus g are in the correspon- 
dence with the divisions of the set of the zeroes of the polynomial P2g+2 on two 
disjoint nonempty subsets with even number of elements. We will consider those 
coverings which correspond to the divisions on subsets with equal number of ele- 
ments and wc can call them equallysplitted, since the Prym variety splits then as a 
sum of two varieties of equal dimension. 

Now, let us consider with the fixed operator A from (5) the whole hierarchy of 
systems defined by the Lax equations 

where 

L(f )(A) = A^S + A^-^Mi + • • • + Mjv 

is a polynomial in A of degree N > 2, and the matrix B is proportional to the 
matrix x- B = dx- From the Lax equations we get the system 

OMn 



dt 

dMk 



+[x,Mfe+i] = [Mfe,f2], 
[x,Mi] = [B,n]. 

Generalizing the situation from the Section 4, we see that the spectral curve Fjv 
is a singular curve of the form 

Pjv(A,m) = + Pjv(A)/x2 + [Qj^{\)f = 0, 

where the polynomials Pn,Qn have degree deg Pjv = degQjv = 2A/'. So, its nor- 
malization is a double covering over the hyperelliptic curve 

i4 = ^ - Ql{X) 

of genus qn = — 1. This covering corresponds to the division of the set of 
zeroes on subsets of zeroes of the polynomials PAr/2 — Qn and PAr/2 + Qn. This 
is an equallysplitted covering under the assumption |xi2| =/= IxsaI we fixed at the 
beginning. It is easy to see that all equallysplitted coverings can be realized in such 
a way. So we have 
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Theorem 3. The Lagrange bitop hierarchy realizes all equally splited coverings over 
the hyperelliptic curves of genus greater than two. 
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